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We construct Z3 vortex solutions in a model in which 5(7(3) is spontaneously broken to Z3. The model is truncated to one 
in which there are only two dimensionless free parameters and the interaction of vortices within this restricted set of models 
is studied numerically. We find that there is a curve in the two dimensional space of parameters for which the energy of two 
asymptotically separated vortices equals the energy of the vortices at vanishing separation. This suggests that the inter-vortex 
potential for Z3 strings might be flat for these couplings, much like the case of (7(1) strings in the Bogomolnyi limit. However, we 
argue that the intervortex potential is attractive at short distances and repulsive at large separations leading to the possibility 
of unstable bound states of Z3 vortices. 



I. INTRODUCTION 

Vortex solutions are well-studied in a wide variety of 
condensed matter systems. In superconductors, vortex 
solutions have been recognized since Abrikosov's seminal 
work H , while string solutions in relativistic field theory 
were first found by Nielsen and Olesen Q. The inter- 
action of vortices has also been a subject of continual 
investigation, starting from the work of Abrikisov who 
conjectured a vortex lattice due to the repulsive force 
between vortices. In relativistic models, early work on 
the interaction of vortices was carried out by Jacobs and 
Rebbi || in which they noted a transition from attractive 
to repulsive interaction as a certain parameter was var- 
ied. Furthermore, at a critical point in parameter space, 
the inter-vortex potential was found to be flat and this 
is the so called Bogomolnyi limit [Q. 

The investigations thus far have mostly considered the 
interaction of C/(l) vortices - the kind that commonly oc- 
cur in superconductors. However, there is a much wider 
variety of vortices occurring in other condensed matter 
systems and there is a possibility that these may also ex- 
ist in particle physics and cosmology. In particular there 
is a class of vortices called "Zn vortices" in which N 
vortices are topologically equivalent to the vacuum. The 
simplest of these is the Z2 (global) vortex that exists in 
nematic liquid crystals. Z4 vortices can be found in the 
A-phase of He 3 . To our knowledge, Z3 vortices have not 
yet been observed but it is possible that these may be 
relevant to confinement in QCD. This is apparent in the 
dual standard model picture that one of us has proposed 
|^| and is indicated by ongoing work on supersymmetric 
dualities ||. 

In this paper we study Z% vortices and their interac- 
tion. The symmetry breaking pattern we consider is 

SU(3) -> Z 3 (1) 

which can be accomplished by the vacuum expectation 
value (VEV) of three adjoint scalar fields. These details 
are provided in Sec. O. The vacuum manifold of the 



model is 

s _ SU(3) 
Z 3 

and since, 

7Tl(£) = Z 3 , 

the model admits Z 3 strings which we explicitly construct 
in Se c. |lH] The interaction of these strings is studied in 
Sec. pyrby comparing the energy of infinitely separated 
(very distant) strings with the energy of the strings at 
vanishing separation. 

An interesting result that we obtain is that there is a 
surface in parameter space such that the energy of two 
vortices at infinite and at vanishing separation are equal. 
This raises the possibility that perhaps the inter-vortex 
potential is flat for these values of the parameters, much 
like the Bogomolnyi case for U(l) strings. However, we 
argue that the Z 3 inter-vortex potential is not flat but 
has a maximum at some finite vortex separation. This 
then indicates that there must exist an unstable, static, 
bound state of two separated Z3 vortices. 

In an early paper, de Vega and Schaposhnik 0] inves- 
tigated the properties of Z3 strings. The group theoretic 
formalism developed there is very general and can be 
used to construct Zn strings for arbitrary N. The study 
of the properties of Z3 strings was, however, restricted to 
a choice of parameters where the Z 3 strings are effectively 
identical to U(l) strings. Furthermore, for this choice of 
parameters, the desired symmetry breaking pattern (|l|) 
is not uniquely picked out by the potential. That is, the 
Z% symmetric vacuum is degenerate with vacuua hav- 
ing other symmetries and is no longer the unique ground 
state of the model. We explain these comments in more 
detail in Sec. [v|. 

II. MODEL 

The construction of a model exhibiting ST/ (3) — > Z3 
involves the following two steps. First, identify the nec- 
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essary ingredients of the model. In particular, determine 
the scalar field content of the model. Second, construct 
the most general scalar field potential and determine the 
range of parameters which lead to the desired symmetry 
breaking. We treat these steps in the next two subsec- 
tions. 

A. Ingredients 

A scalar field, <&, in the adjoint representation of SU (3) 
can be written as the 3x3 Hermitian matrix 

8 

$ = ]T $ a A a 

a=l 

where, A a are the Gell-Mann matrices l| and <E> a are 
8 real scalar fields. Under the action of g G SU(3), $ 
transforms as: 

$ — * $' = g^g' 1 , 

where, g may be written as: exp(ia a A a ) for any set of 
a a . 

Now the centeiQ of SU(3) is Z 3 and the elements of Z 3 
are of the form: 

e i27rn/3 1 ^ n _ mr ,gg er j 

where 1 is the identity matrix. Hence, $ is left invariant 
under transformations belonging to the center. And so 
the VEVs of any number of adjoint scalar fields cannot 
break the Z 3 center. Then, one way to achieve ((!]) would 
be to give VEVs to as many adjoint fields as necessary 
to break the SU(3) maximally, that is, only leaving the 
center unbroken. Indeed, one can check that two adjoints 
$1 and $2 are sufficient because if, 

$1 = Ai , $ 2 = A 4 (2) 

then the only group elements that commute with both 
$1 and $2 are the ones proportional to 1, that is, the 
elements of the Z 3 center. For completeness, in Table | 
we show the various possible VEVs for ($1, $2) and the 
resulting symmetry breaking pattern. 

In the next subsection, we will consider an SU (3) in- 
variant potential for two adjoint scalar fields. As the full 
potential has a lot of parameters, we restrict our atten- 
tion to a certain region of parameter space. Within this 
truncated model, we will find that it is not possible to 
construct a potential that will lead to the VEVs in (||) at 
a unique, global minimum. Then it will be necessary to 



introduce a third adjoint field $3. In Table || we show 
the symmetry breaking patterns for different directions 

of (<!>!, $ 2 , $3). 

In what follows, we will construct a potential that will 
have a global minimum when ($1, $2, $3) acquire VEVs 
in the following direction: 

(A4, A 6 , A 8 ) . (3) 

(SU(3) rotations of these VEVs will yield the manifold 
of global minima.) As shown in Table II, the residual 
symmetry group will be Z 3 in this case. 
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TABLE I. SU(3) breaking with two adjoint Higgs fields. 
The rows and columns label the direction of the VEVs of 
each of the two fields while the table entry gives the resid- 
ual symmetry. For convenience of notation we have defined: 
(71(71 = (7(1) x U(l), and SU2U1 = 5(7(2) x (7(1) 



U(l) 


(Xi, X2, X 3 ), (Xi, X2, Xs), (Ai,A3,A8), 
(A2,A3,As), (A3,A4,As), (A3,A5,As), 
(A3,A6,As), (A3,A7,As), (A4,A5,As), 
(A6, A7, Ag), 


z 3 


all other (Xi,Xj,Xk) with distinct i,j,k 



TABLE II. 5(7(3) symmetry breaking with three distinct 
adjoint scalar fields. The first column shows the residual sym- 
metry group if the fields get VEVs in the directions shown in 
the second column. 



*The center of a group consists of elements that commute 
with all other elements. 
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B. Construction of the potential 



and, make the further choice of parameters: 



The SU(3) invariant potential for three adjoint fields 
can be written as follows: 



V({$i}) = V 1 ({$ 1 }) + V 2 ({$t}) 



(4) 



where 



Vi({*i}) = ]T[-m 2 (Tr<i> 2 ) + ai (Tr^) 2 + 6 ; Tr$f] (5) 
1=1 

and 

3 

V2({$j» = 5>,Tr($ m $„) + d / (Tr($ m $„)) 2 (6) 



+e,Tr($ m <i> n ) 2 + fiTr(<S> 2 m <S>l) 
+ ffi Tr($ 2 $ m $„) + ^Tr($ m $ 2 $„)] 

with, l,m,n taking cyclic values over 1,2,3. (In writing 
the potential we have omitted cubic terms for simplicity.) 
The full Lagrangian can now be written: 

L = J2 Tr(|^$/| 2 ) - ^TriG^On - V({$,}) (7) 



where, 



D^i = + giA^^i] , 



is the matrix-valued gauge field, 

Gfiv = G^„A a = d^Au — d v A^ + g[A^, A v \ . 

The general potential (Q) has 27 parameters and is far 
too complicated for us to handle. We will assume cer- 
tain relationships among the parameters to enable us to 
proceed further. These are: 



ai = a , hi = b , m; 



/?? 



and 



Q = = di = gi = hi . 



Now to check if two scalar fields would have been suf- 
ficient for the symmetry breaking ([!]), we eliminate all 
terms containing $3 in ([!]), restrict our attention to the 
truncated region in parameter space, then feed in the 
various possible VEVs for $1 and $2 from Table [j]. We 
find that the VEVs leading to a Z 3 residual symmetry 
give a higher (or equal) energy than the other symmetry 
breaking patterns. So the Z% vacuum cannot be a unique 
global minimum. Hence, it is necessary for us to include 
the third scalar field $ 3 . 

With three scalar fields, we introduce 

A = 16a + 8b , 



fl = ei = e-2 



-e3 



where, e is a new free parameter. The minus sign in 
front of e 3 is a crucial feature that ensures that the global 
minimum of the potential has Z 3 symmetry 

The potential for this restricted set of parameters is 



3 

£ 

1=1 



-(Tr$f) + a(Tr$f) 2 + oTr$f] 



+e [ Tr($ 2 $ 3 ) 2 + Tr($ 3 $!) 2 - Tr($!$ 2 ) 2 
+ Tr($ 2 $ 2 ) + Tr($ 2 $ 2 ) + Tr($ 2 $ 2 )] . (8) 

For the potential to have a global minimum at finite 
VEVs of the fields, we need 

A > . 

The requirement that the VEVs of the fields be non- 
vanishing gives the constraint: 

e 3 
A<2- 

Within this parameter range, we have inserted all possi- 
ble choices of directions of $/ (I = 1, 2, 3) in the potential 
and find that it has a global minimum when the residual 
symmetry group is Z3, provided 

e > . 

Further, the VEVs yielding the global minimum are in 
the (A4,A6,Ag) directions (and SU(3) rotations of these 
directions): 



$1 = ??iA 4 
$2 = ??2A 6 

$3 = %As 



(9) 



where 



m=m = m\ 



m = m \ 



A- 2e/3 



A 2 + 2eA - 8e 2 /9 



A + 2e/3 



A 2 + 2eA - 8e 2 /9 



(10) 



(11) 



Let us now define: 



-Tr($ 



Then, for a vortex, vi will vary in space and the relevant 
potential is: 



V({vi}) = J2 



A 



"'7 



+ ^(vi + v 2 2 )v 2 3 + ev 2 lV 2 2 . (12) 
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III. Z z STRING ANSATZ AND SOLUTION 

Having truncated the full model to one which is simple 
enough to analyze, we now write down the ansatz for Z% 
strings, insert it into the field equations and then find the 
string solutions. 

To write a string ansatz, we must first specify a closed 
path on the vacuum manifold, P(9), parametrized by 9 S 
[0, 2ir], which is incontractable. This is given by: 

P{9) = e inX sO/V3 

This path is incontractable since P(2tt) is a non-trivial 
element of the discrete residual group Z 3 for n — 1,2. 

We now identify 9 with the spatial polar coordinate. 
Then the scalar field ansatz is: 

$i(r -kx>,6) = P(9)^i(r oo, = O)P(0) . 

With $i(6 = 0) given by eq. @, this leads to: 

<fr 1 (r, 9) = v 1 (r)(cosn9 A 4 + sinn9 A 5 ) (13) 
$2(r,9) — V2(r)(cosn9 Xq + sinn9 A7) 
®3{r,6) = v 3 {r)X 8 
A s n a ( r ) 

where we have included scalar field profile functions vi (r) 
and given the gauge field ansatz with its profile function 
a(r). All other components of the gauge field are taken 
to vanish. 

We now insert this ansatz into the field equations to 
get: 



1 



(1 - a) 2 vi - Auf + m 2 vi 
1 

— 1 
3 



-M4 



-v' 2 J C 1 _ a f V 1 ~ AlJ 2 + m2 «2 



-2e(vl + 3«3)«2 





A«3 + m A v 3 — -e(«x + v 2) v 3 = 
- -a' + 3g 2 {v 2 + vl){\ - a) = 



' • ..3 , „,2„, ^ , 1 2 



a 



where a prime denotes differentiation with respect to r. 

The equations for v± and ^2 are identical and so are 
the asymptotic boundary conditions as is seen from (10). 
Therefore we set 

vi(r) = v 2 (r) . 

It is now convenient to define rescaled coordinates, pa- 
rameters and fields as follows: 



'6 

' —mr 
A 



(14) 



A 



3 5 2 



777. Til 
Vl=V 2 = -^f( x ) > v 3 = ~Jj h ( X } 

Primes will now denote differentiation with respect to the 
rescaled coordinate x. 

The rescaled equations are: 



r+^-i(i-«) 2 /-^-(/ 2 -D/ 



A 2 
2 



P(f + -5-)/ = (15) 



a; 2 



a "__ + /2(i_ a ) = 



(16) 



(17) 



The boundary conditions on the functions /, h and 
a follow by requiring single- valuedness and regularity of 
the fields at the origin, 



«(0) = /(0) = h'(0) = 



(18) 



At infinity, the fields should go to their vacuum expecta- 
tion values: 

a(x — > 00) = 1 , /(x — > cxd) = _F 1 M 2 - ~ * 00 ) = ^0 

(19) 

with, 



_ 1-2,7/3 

y 1 + 2,7-8^/9 ' ^° 



' 1 + 2cr/3 
l + 2cr-8a 2 /9 ' 

(20) 



where 

CT ~ A 2 ' 

If we set /3 = 0, the ft. equation is solved by h = 1 and 
the / and a equations are exactly the Nielson-Olesen 
equations for the Abelian-Higgs vortex. These vortices 
have been studied extensively and this is also the case 
discussed by de Vega and Schaposhnik |Q. The interac- 
tion of Abelian-Higgs vortices is characterized by the sin- 
gle ratio of length scales entering the problem - namely, 
the ratio of the (single) scalar and vector masses. With 
(3 7^ 0, however, the picture is more complicated since 
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we have three length scales corresponding to the masses 
of the two independent scalar fields in our truncated 
model, and the gauge fields. Therefore there are two 
independent dimensionless ratios we can construct, and 
a correspondingly richer structure to the interaction of 
Z3 vorticeal. 

We have solved the equations of motion (|l5])-(|l7|) by 
a numerical shooting routine. In Fig. [I] we show the 
behaviour of the /, h and v fields for a particular choice 
of parameters and for n = 1. 

FIG. 1. The profile functions /, h and a for the unit wind- 
ing Z-i vortex for A = 1.13 and /3 = 0.42. 
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The equations of motion for the n — 2 vortex can also 
be solved. So the n = 2 vortex is indeed a solution, 
though it will be unstable to decay into the topologically 
equivalent n = — 1 solution which has lower energy. This 
instability of the n — 2 vortex is not relevant for us since 
all that we are interested in is the energy of two overlap- 
ping n=l vortices which is the same as the energy of 
the n = 2 vortex solution. 



IV. ENERGY FUNCTIONAL AND NUMERICAL 
EVALUATION 

The set of field equations (H)-© can also be ob- 
tained by extremizing the energy functional 

E = Jd 3 x [\G%G aij + i(A*?)(A*f) + v({^})} 

(21) 

(The index a = 1, 8 is the group index and I = 1, 2, 3 
labels the different adjoint fields.) Here we will obtain 



^In the full model with 27 parameters, one might expect the 
interaction to be even more complex. However, our under- 
standing of the results presented in Sec. IVD indicate that the 
picture is likely to be quite simple even in the full model. 



the extremum values of the energy for the n = 1, 2 topo- 
logical configurations directly, that is, without solving 
the field equations of motion. This is the technique used 
by Jacobs and Rebbi H to study the interaction of U(l) 
vortices, and we will employ it to study the interaction 
of Z3 vortices. 
Let us first define 



a = X\ 1 <J , v = Ai / ^— 

1 3 ' V l + 2cr -8cr 2 /9 



(22) 



Then the solution to eq. (|l5|)-([l~7|) may be written as 



f = F 1 + £ 



3- 



3=0 



3- 



(23) 



(24) 



(25) 



In this form, the functions automatically satisfy the de- 
sired boundary at infinity (eq. ([l~9|)). The boundary con- 
ditions at the origin (eq. (|l8|)) require 



a 



-1 , fo 



-1 , h\ = lyfiQ 



for both 7i = 1 and n — 2 vortices. In addition, regularity 
of the gauge fields at the origin requires 



a'(0) = 



and, 



/'(0) = 0, for n = 2 . 
These conditions give 

ai = —1 

and 

fi = -[I, for n = 2 . 

In this scheme, since we are numerically evaluating the 
energy functional, it is necessary to ensure that the po- 
tential vanishes in the true vacuum. This requires that 
we shift the potential in ( |l~2| ) by a constant Vq: 

v = -(2F 4 + H*) + 2(2F 2 + H 2 a ) - + 4aF 4 . 

The potential can now be written in terms of the fields 
/ and h: 
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\2 r 



V(f,h) = - (2/ 4 + h 4 ) - 2(2 f + h 2 ) 



+ -afh 2 + 4af 4 + vo 



In terms of the fields /, h and a, the energy is: 

27TTO 



E = 



6Ag 



J dz dx x £[f, h, a] 



where z is the rescaled (as in eq. ) dimensionless 
coordinate along the string, and 

£[f, h, a] = f' 2 + \h n + n 2 f + 4(1 - «)7 2 
2 \ x J x 

+ V(f,h) . 

We now have to evaluate the energy functional in terms 
of the various coefficients /j, Tij and Oj (infinite in num- 
ber), for different choices of {A,/3} and winding number 
n. On inserting the expansions for the fields in the en- 
ergy functional, we end up with integrals that can be 
evaluated in terms of Gamma functions. These form the 
coefficients of the quartic order polynomials in the fi, 
hi and a*. This gives E = E({fi}, {hi}, {at}) where 
i = 0, ...,oo. We then truncate the expansion by re- 
taining only the first 7 terms in each of the expansions 
(i = 0, ...,6), making a total of 21 parameters that need 
to be varied to minimize E. We finally find the global 
minimum of E with respect to the variation of the coef- 
ficients for n = 1, 2 and A G (0, 10) and (3 G (0, 6). 

Sample results for the dependence of the energy on the 
/3 parameter are shown in Fig. ^. 

FIG. 2. The dependence of the energy of n = 1 and n — 2 
vortices on the parameter f3 for A = 2.0. Also shown is twice 
the energy of the n = 1 vortex. 
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V. CRITICAL COUPLINGS AND 
INTER- VORTEX FORCES 

In the = limit, the equation of motion ( |l6| ) is 
simply solved by 

h(x;f3 = 0) = 1 

and the remaining equations (|l5|) and ( |l7|) are identical 
to the U(l) equations. So, in this case, the structure and 
interaction of Z 3 vortices is identical to those of 17(1) vor- 
tices. In particular, for A = 1, the inter- vortex potential 
vanishes [Q. If, however, we now return to the potential 
(||) and set e = (equivalent to ft = 0), we find that the 
potential does not uniquely pick out the required direc- 
tions for the symmetry breaking (Q). Indeed, the VEVs 
of the three different fields can point in the same direc- 
tion, say in the As direction, leading to an SU(2) x U(l) 
residual symmetry. So to pick out the desired symmetry 
breaking we must necessarily consider (3 =/= 0. 

As one can see in Fig. p|, there are points in parameter 
space where 



E(/3, A; n = 2) = 2 x E(J3, A; n = 1) 



(26) 



In this case, the energy of two infinitely separated vor- 
tices is equal to the energy of two overlapping vortices. 
We shall call the parameters for which (|2^) holds, to be 
"critical" . In Fig. || we plot the critical curve in (/3, A) 
space. 

FIG. 3. The curve in parameter space for which the energy 
of two inifintely separated vortices have the same energy as 
two vortices at zero separation, that is, the n — 2 vortex. The 
region in which the n = 2 vortex is more energetic than two 
n = 1 vortices is the "repulsive" region while that in which 
the n — 2 vortex is less energetic is the "attractive" region. 




attractive region 



beta 
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The question we now address is whether the critically 
coupled vortices have a mutual repulsion or attraction at 
intermediate separations. In the critically coupled U(l) 
case (ft — 0, A = 1), the intervortex potential is flat and 
hence the intervortex forces vanish at all separations. It is 
useful to think of this in terms of the repulsive gauge field 
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and attractive scalar field interactions. The exchange of 
spin one gauge particles between identical vortices leads 
to a repulsive force while the exchange of spin zero par- 
ticles leads to an attractive force. So the gauge field 
repulsion is balanced by the scalar field attraction in the 
critical U(l) case. 

In the case of the Z3 string, the crucial observation is 
that the curve of critical couplings lies in the region A > 1 
and small a — (3/\ 2 - Hence, along the critical curve, the 
gauge field mass is smaller than either of the scalar field 
masses /1 and v given in ( p2| ) . So the gauge field interac- 
tion is of longer range than the scalar field interaction. 
Therefore, if we bring in two infinitely separated vor- 
tices, they will first experience the repulsive gauge field 
interaction. When they come in closer, the attractive 
interaction due to the scalar fields will turn on. Since 
there are two scalar fields in our model, the attraction is 
stronger than in the U(l) case and is more effective in 
cancelling out the gauge field repulsion. So the intervor- 
tex potential is expected to turn over, as schematically 
depicted in Fig. For the critically coupled case, the 
turn over is such that the energy at zero vortex separation 
equals that at infinite vortex separation. The presence 
of a turning point in the intervortex potential at some 
vortex separation s* means that two vortices separated 
by this distance can be in relative equilibrium. In the 
present case, this is an unstable equilibrium because the 
potential is a maximum. 

FIG. 4. A schematic depiction of the expected intervortex 
potential (U) as a function of vortex separation (s) in the 
critically coupled case where the gauge field mass is less than 
the scalar field masses. Being lighter, the gauge field provides 
a longer range interaction than the scalar fields and leads to 
a repulsive force between vortices in the asymptotic region. 
In a U(l) model, the single scalar field is unable to overcome 
the repulsive potential at short distances. However, in the Zz 
string case, the two scalar fields successfully turn the potential 
over at short distances leading to a maximum at s = s*. 

U 



feature can be understood by realizing that the presence 
of two scalar fields instead of one implies an enhanced 
attractive force between vortices at short distances. So 
if one imagines starting out with repulsive U(l) vortices, 
and adding a second scalar field as is present in the Z3 
vortices, this can turn over the intervortex potential at 
small separations and make the energy at s = the same 
as that at s — 00. But since one needs to start out 
with repulsive U(l) vortices, this means that critically 
coupled Z3 vortices should lie in the A > 1 region. (If 
one started out with attractive (A < 1) U(l) vortices, 
the addition of a second scalar field would simply make 
the vortices even more attractive at short distances and 
eq. ( ^6|) could never be satisfied.) 

VI. CONCLUSIONS 

We have constructed field theory solutions for Z3 
strings and have studied their interaction within a range 
of model parameters. The solution for the structure of 
the vortex is shown in Fig. [I] while the dependence of 
the energy on the new parameter in the model, that is, 
the parameter not present in the U(l) case, is shown in 
Fig. 0. Infinitely separated Z3 strings have the same 
energy as strings at zero separation along a curve in our 
two dimensional space of parameters (Fig. ||). However, 
we have given general arguments to show that, unlike the 
U(l) case, the intervortex potential is not trivial for these 
critically coupled Z3 strings. In fact, the intervortex po- 
tential is expected to have a maximum value at some 
non- vanishing vortex separation (Fig. |4|). This suggests 
that an unstable bound state of two Z3 vortices should 
exist. 

Although we have worked in detail within a specific 
range of parameters, we have understood the intervortex 
forces based on the number of scalar and vector fields 
present in the model and their masses. This reasoning 
(described in Sec. |v|) is expected to be valid quite gener- 
ally and should apply to the full range of parameters in 
this model as well as to other models. 
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Note that this argument is based on the fact that the 
curve of critical couplings lies in the A > 1 region. This 
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